A large normative connectome for exploring the tractographic correlates of focal brain interventions.

Diffusion-weighted MRI (dMRI) is a widely used neuroimaging modality that permits the in vivo exploration of white matter connections in the human brain. Normative structural connectomics - the application of large-scale, group-derived dMRI datasets to out-of-sample cohorts - have increasingly been leveraged to study the network correlates of focal brain interventions, insults, and other regions-of-interest (ROIs). Here, we provide a normative, whole-brain connectome in MNI space that enables researchers to interrogate fiber streamlines that are likely perturbed by given ROIs, even in the absence of subject-specific dMRI data. Assembled from multi-shell dMRI data of 985 healthy Human Connectome Project subjects using generalized Q-sampling imaging and multispectral normalization techniques, this connectome comprises ~12 million unique streamlines, the largest to date. It has already been utilized in at least 18 peer-reviewed publications, most frequently in the context of neuromodulatory interventions like deep brain stimulation and focused ultrasound. Now publicly available, this connectome will constitute a useful tool for understanding the wider impact of focal brain perturbations on white matter architecture going forward.

Open access EEG dataset of repeated measurements from a single subject for microstate analysis.

Electroencephalography (EEG) microstate analysis is a neuroimaging analytical method that has received considerable attention in recent years and is widely used for analysing EEG signals. EEG is easily influenced by internal and external factors, which can affect the repeatability and stability of EEG microstate analysis. However, there have been few reports and publicly available datasets on the repeatability of EEG microstate analysis. In the current study, a 39-year-old healthy male underwent a total of 60 simultaneous electroencephalography and electrocardiogram measurements over a period of three months. After the EEG recording was completed, magnetic resonance imaging (MRI) was also conducted. To date, this EEG dataset has the highest number of repeated measurements for one individual. The dataset can be used to assess the stability and repeatability of EEG microstates and other analytical methods, to decode resting EEG states among subjects with open eyes, and to explore the stability and repeatability of cortical spatiotemporal dynamics through source analysis with individual MRI.

Multiomics analysis to explore blood metabolite biomarkers in an Alzheimer's Disease Neuroimaging Initiative cohort.

Alzheimer's disease (AD) is a neurodegenerative disease that commonly causes dementia. Identifying biomarkers for the early detection of AD is an emerging need, as brain dysfunction begins two decades before the onset of clinical symptoms. To this end, we reanalyzed untargeted metabolomic mass spectrometry data from 905 patients enrolled in the AD Neuroimaging Initiative (ADNI) cohort using MS-DIAL, with 1,304,633 spectra of 39,108 unique biomolecules. Metabolic profiles of 93 hydrophilic metabolites were determined. Additionally, we integrated targeted lipidomic data (4873 samples from 1524 patients) to explore candidate biomarkers for predicting progressive mild cognitive impairment (pMCI) in patients diagnosed with AD within two years using the baseline metabolome. Patients with lower ergothioneine levels had a 12% higher rate of AD progression with the significance of P = 0.012 (Wald test). Furthermore, an increase in ganglioside (GM3) and decrease in plasmalogen lipids, many of which are associated with apolipoprotein E polymorphism, were confirmed in AD patients, and the higher levels of lysophosphatidylcholine (18:1) and GM3 d18:1/20:0 showed 19% and 17% higher rates of AD progression, respectively (Wald test: P = 3.9 × 10-8 and 4.3 × 10-7). Palmitoleamide, oleamide, diacylglycerols, and ether lipids were also identified as significantly altered metabolites at baseline in patients with pMCI. The integrated analysis of metabolites and genomics data showed that combining information on metabolites and genotypes enhances the predictive performance of AD progression, suggesting that metabolomics is essential to complement genomic data. In conclusion, the reanalysis of multiomics data provides new insights to detect early development of AD pathology and to partially understand metabolic changes in age-related onset of AD.

Neuroimaging alterations and relapse in early-stage psychosis.

BACKGROUND: Recent reports have indicated that symptom exacerbation after a period of improvement, referred to as relapse, in early-stage psychosis could result in brain changes and poor disease outcomes. We hypothesized that substantial neuroimaging alterations may exist among patients who experience relapse in early-stage psychosis. METHODS: We studied patients with psychosis within 2 years after the first psychotic event and healthy controls. We divided patients into 2 groups, namely those who did not experience relapse between disease onset and the magnetic resonance imaging (MRI) scan (no-relapse group) and those who did experience relapse between these 2 timings (relapse group). We analyzed 3003 functional connectivity estimates between 78 regions of interest (ROIs) derived from resting-state functional MRI data by adjusting for demographic and clinical confounding factors. RESULTS: We studied 85 patients, incuding 54 in the relapse group and 31 in the no-relapse group, along with 94 healthy controls. We observed significant differences in 47 functional connectivity estimates between the relapse and control groups after multiple comparison corrections, whereas no differences were found between the no-relapse and control groups. Most of these pathological signatures (64%) involved the thalamus. The Jonckheere-Terpstra test indicated that all 47 functional connectivity changes had a significant cross-group progression from controls to patients in the no-relapse group to patients in the relapse group. LIMITATIONS: Longitudinal studies are needed to further validate the involvement and pathological importance of the thalamus in relapse. CONCLUSION: We observed pathological differences in neuronal connectivity associated with relapse in early-stage psychosis, which are more specifically associated with the thalamus. Our study implies the importance of considering neurobiological mechanisms associated with relapse in the trajectory of psychotic disorders.

Computational limits to the legibility of the imaged human brain.

Our knowledge of the organisation of the human brain at the population-level is yet to translate into power to predict functional differences at the individual-level, limiting clinical applications and casting doubt on the generalisability of inferred mechanisms. It remains unknown whether the difficulty arises from the absence of individuating biological patterns within the brain, or from limited power to access them with the models and compute at our disposal. Here we comprehensively investigate the resolvability of such patterns with data and compute at unprecedented scale. Across 23 810 unique participants from UK Biobank, we systematically evaluate the predictability of 25 individual biological characteristics, from all available combinations of structural and functional neuroimaging data. Over 4526 GPU*hours of computation, we train, optimize, and evaluate out-of-sample 700 individual predictive models, including fully-connected feed-forward neural networks of demographic, psychological, serological, chronic disease, and functional connectivity characteristics, and both uni- and multi-modal 3D convolutional neural network models of macro- and micro-structural brain imaging. We find a marked discrepancy between the high predictability of sex (balanced accuracy 99.7%), age (mean absolute error 2.048 years, R2 0.859), and weight (mean absolute error 2.609Kg, R2 0.625), for which we set new state-of-the-art performance, and the surprisingly low predictability of other characteristics. Neither structural nor functional imaging predicted an individual's psychology better than the coincidence of common chronic disease (p < 0.05). Serology predicted chronic disease (p < 0.05) and was best predicted by it (p < 0.001), followed by structural neuroimaging (p < 0.05). Our findings suggest either more informative imaging or more powerful models will be needed to decipher individual level characteristics from the human brain. We make our models and code openly available.

Sex, gender diversity, and brain structure in early adolescence.

There remains little consensus about the relationship between sex and brain structure, particularly in early adolescence. Moreover, few pediatric neuroimaging studies have analyzed both sex and gender as variables of interest-many of which included small sample sizes and relied on binary definitions of gender. The current study examined gender diversity with a continuous felt-gender score and categorized sex based on X and Y allele frequency in a large sample of children ages 9-11 years old (N = 7195). Then, a statistical model-building approach was employed to determine whether gender diversity and sex independently or jointly relate to brain morphology, including subcortical volume, cortical thickness, gyrification, and white matter microstructure. Additional sensitivity analyses found that male versus female differences in gyrification and white matter were largely accounted for by total brain volume, rather than sex per se. The model with sex, but not gender diversity, was the best-fitting model in 60.1% of gray matter regions and 61.9% of white matter regions after adjusting for brain volume. The proportion of variance accounted for by sex was negligible to small in all cases. While models including felt-gender explained a greater amount of variance in a few regions, the felt-gender score alone was not a significant predictor on its own for any white or gray matter regions examined. Overall, these findings demonstrate that at ages 9-11 years old, sex accounts for a small proportion of variance in brain structure, while gender diversity is not directly associated with neurostructural diversity.

ANTsX neuroimaging-derived structural phenotypes of UK Biobank.

UK Biobank is a large-scale epidemiological resource for investigating prospective correlations between various lifestyle, environmental, and genetic factors with health and disease progression. In addition to individual subject information obtained through surveys and physical examinations, a comprehensive neuroimaging battery consisting of multiple modalities provides imaging-derived phenotypes (IDPs) that can serve as biomarkers in neuroscience research. In this study, we augment the existing set of UK Biobank neuroimaging structural IDPs, obtained from well-established software libraries such as FSL and FreeSurfer, with related measurements acquired through the Advanced Normalization Tools Ecosystem. This includes previously established cortical and subcortical measurements defined, in part, based on the Desikan-Killiany-Tourville atlas. Also included are morphological measurements from two recent developments: medial temporal lobe parcellation of hippocampal and extra-hippocampal regions in addition to cerebellum parcellation and thickness based on the Schmahmann anatomical labeling. Through predictive modeling, we assess the clinical utility of these IDP measurements, individually and in combination, using commonly studied phenotypic correlates including age, fluid intelligence, numeric memory, and several other sociodemographic variables. The predictive accuracy of these IDP-based models, in terms of root-mean-squared-error or area-under-the-curve for continuous and categorical variables, respectively, provides comparative insights between software libraries as well as potential clinical interpretability. Results demonstrate varied performance between package-based IDP sets and their combination, emphasizing the need for careful consideration in their selection and utilization.

Joint transformer architecture in brain 3D MRI classification: its application in Alzheimer's disease classification.

Alzheimer's disease (AD), a neurodegenerative disease that mostly affects the elderly, slowly impairs memory, cognition, and daily tasks. AD has long been one of the most debilitating chronic neurological disorders, affecting mostly people over 65. In this study, we investigated the use of Vision Transformer (ViT) for Magnetic Resonance Image processing in the context of AD diagnosis. ViT was utilized to extract features from MRIs, map them to a feature sequence, perform sequence modeling to maintain interdependencies, and classify features using a time series transformer. The proposed model was evaluated using ADNI T1-weighted MRIs for binary and multiclass classification. Two data collections, Complete 1Yr 1.5T and Complete 3Yr 3T, from the ADNI database were used for training and testing. A random split approach was used, allocating 60% for training and 20% for testing and validation, resulting in sample sizes of (211, 70, 70) and (1378, 458, 458), respectively. The performance of our proposed model was compared to various deep learning models, including CNN with BiL-STM and ViT with Bi-LSTM. The suggested technique diagnoses AD with high accuracy (99.048% for binary and 99.014% for multiclass classification), precision, recall, and F-score. Our proposed method offers researchers an approach to more efficient early clinical diagnosis and interventions.

Beyond visual integration: sensitivity of the temporal-parietal junction for objects, places, and faces.

One important role of the TPJ is the contribution to perception of the global gist in hierarchically organized stimuli where individual elements create a global visual percept. However, the link between clinical findings in simultanagnosia and neuroimaging in healthy subjects is missing for real-world global stimuli, like visual scenes. It is well-known that hierarchical, global stimuli activate TPJ regions and that simultanagnosia patients show deficits during the recognition of hierarchical stimuli and real-world visual scenes. However, the role of the TPJ in real-world scene processing is entirely unexplored. In the present study, we first localized TPJ regions significantly responding to the global gist of hierarchical stimuli and then investigated the responses to visual scenes, as well as single objects and faces as control stimuli. All three stimulus classes evoked significantly positive univariate responses in the previously localized TPJ regions. In a multivariate analysis, we were able to demonstrate that voxel patterns of the TPJ were classified significantly above chance level for all three stimulus classes. These results demonstrate a significant involvement of the TPJ in processing of complex visual stimuli that is not restricted to visual scenes and that the TPJ is sensitive to different classes of visual stimuli with a specific signature of neuronal activations.

A bilateral filtering-based image enhancement for Alzheimer disease classification using CNN.

This study aims to develop an optimally performing convolutional neural network to classify Alzheimer's disease into mild cognitive impairment, normal controls, or Alzheimer's disease classes using a magnetic resonance imaging dataset. To achieve this, we focused the study on addressing the challenge of image noise, which impacts the performance of deep learning models. The study introduced a scheme for enhancing images to improve the quality of the datasets. Specifically, an image enhancement algorithm based on histogram equalization and bilateral filtering techniques was deployed to reduce noise and enhance the quality of the images. Subsequently, a convolutional neural network model comprising four convolutional layers and two hidden layers was devised for classifying Alzheimer's disease into three (3) distinct categories, namely mild cognitive impairment, Alzheimer's disease, and normal controls. The model was trained and evaluated using a 10-fold cross-validation sampling approach with a learning rate of 0.001 and 200 training epochs at each instance. The proposed model yielded notable results, such as an accuracy of 93.45% and an area under the curve value of 0.99 when trained on the three classes. The model further showed superior results on binary classification compared with existing methods. The model recorded 94.39%, 94.92%, and 95.62% accuracies for Alzheimer's disease versus normal controls, Alzheimer's disease versus mild cognitive impairment, and mild cognitive impairment versus normal controls classes, respectively.

Mapping of Alzheimer's disease related data elements and the NIH Common Data Elements.

BACKGROUND: Alzheimer's Disease (AD) is a devastating disease that destroys memory and other cognitive functions. There has been an increasing research effort to prevent and treat AD. In the US, two major data sharing resources for AD research are the National Alzheimer's Coordinating Center (NACC) and the Alzheimer's Disease Neuroimaging Initiative (ADNI); Additionally, the National Institutes of Health (NIH) Common Data Elements (CDE) Repository has been developed to facilitate data sharing and improve the interoperability among data sets in various disease research areas. METHOD: To better understand how AD-related data elements in these resources are interoperable with each other, we leverage different representation models to map data elements from different resources: NACC to ADNI, NACC to NIH CDE, and ADNI to NIH CDE. We explore bag-of-words based and word embeddings based models (Word2Vec and BioWordVec) to perform the data element mappings in these resources. RESULTS: The data dictionaries downloaded on November 23, 2021 contain 1,195 data elements in NACC, 13,918 in ADNI, and 27,213 in NIH CDE Repository. Data element preprocessing reduced the numbers of NACC and ADNI data elements for mapping to 1,099 and 7,584 respectively. Manual evaluation of the mapping results showed that the bag-of-words based approach achieved the best precision, while the BioWordVec based approach attained the best recall. In total, the three approaches mapped 175 out of 1,099 (15.92%) NACC data elements to ADNI; 107 out of 1,099 (9.74%) NACC data elements to NIH CDE; and 171 out of 7,584 (2.25%) ADNI data elements to NIH CDE. CONCLUSIONS: The bag-of-words based and word embeddings based approaches showed promise in mapping AD-related data elements between different resources. Although the mapping approaches need further improvement, our result indicates that there is a critical need to standardize CDEs across these valuable AD research resources in order to maximize the discoveries regarding AD pathophysiology, diagnosis, and treatment that can be gleaned from them.

A Multiscale Closed-Loop Neurotoxicity Model of Alzheimer's Disease Progression Explains Functional Connectivity Alterations.

The accumulation of amyloid-ß (Aß) and hyperphosphorylated-tau (hp-tau) are two classical histopathological biomarkers in Alzheimer's disease (AD). However, their detailed interactions with the electrophysiological changes at the meso- and macroscale are not yet fully understood. We developed a mechanistic multiscale model of AD progression, linking proteinopathy to its effects on neural activity and vice-versa. We integrated a heterodimer model of prion-like protein propagation and a brain network model of Jansen-Rit neural masses derived from human neuroimaging data whose parameters varied due to neurotoxicity. Results showed that changes in inhibition guided the electrophysiological alterations found in AD, and these changes were mainly attributed to Aß effects. Additionally, we found a causal disconnection between cellular hyperactivity and interregional hypersynchrony contrary to previous beliefs. Finally, we demonstrated that early Aß and hp-tau depositions' location determine the spatiotemporal profile of the proteinopathy. The presented model combines the molecular effects of both Aß and hp-tau together with a mechanistic protein propagation model and network effects within a closed-loop model. This holds the potential to enlighten the interplay between AD mechanisms on various scales, aiming to develop and test novel hypotheses on the contribution of different AD-related variables to the disease evolution.

The association between cannabis use and neuroimaging measures in older adults: findings from the UK biobank.

BACKGROUND: Cannabis use has increased in recent years. However, the long-term implications of cannabis use on brain health remain unknown. We explored the associations of cannabis use with volumetric brain magnetic resonance imaging (MRI) measures in dementia-free older adults. METHODS: This cross-sectional and longitudinal study included dementia-free participants of the UK Biobank aged ≥60 years. Linear regression models were used to evaluate the association of cannabis use and patterns of use with volumetric brain MRI measures. The association between cannabis use and change in brain MRI measures over time was also tested. All models were adjusted for potential confounders. RESULTS: The sample included 19,932 participants (mean age 68 ± 5 years, 48% men), 3,800 (19%) reported lifetime use of cannabis. Cannabis use was associated with smaller total, white, grey and peripheral cortical grey matter volumes (B = -6,690 ± 1,157; P < 0.001, B = -4,396 ± 766; P < 0.001, B = -2,140 ± 690; P = 0.002 and B = -2,451 ± 606; P < 0.001, respectively). Among cannabis users, longer duration of use was associated with smaller total brain, grey and cortical grey matter volumes (B = -7,878 ± 2,396; P = 0.001, B = -5,411 ± 1,430; P < 0.001, B = -5,396 ± 1,254; P < 0.001, respectively), and with increased white matter hyperintensity volume (B = 0.09 ± 0.03; P = 0.008). Additionally, current vs. former users (B = -10,432 ± 4,395; P = 0.020) and frequent versus non-frequent users (B = -2,274 ± 1,125; P = 0.043) had smaller grey and cortical grey matter volumes, respectively. No significant associations were observed between cannabis use and change in brain MRI measures. DISCUSSION: Our findings suggest that cannabis use, particularly longer duration and frequent use, may be related to smaller grey and white matter volumes in older ages, but not to late-life changes in these measures over time.

Lethal mycotic pseudoaneurysm presenting as isolated sixth nerve palsy.

Pseudoaneurysm of the internal carotid artery caused by skull base osteomyelitis (SBO) is a lethal condition seen in immunocompromised patients, predominantly those with diabetes mellitus. Cranial nerve involvement is a common complication and generally indicates a poor prognosis. We report the case of a 62-year-old diabetic patient who presented with isolated sixth cranial nerve palsy. She had uncontrolled blood sugar levels and high erythrocyte sedimentation rate, and she suffered from pyelonephritis. Neuroimaging detected SBO with multiple secondary mycotic pseudoaneurysms prominent at the petrocavernous junction. Ischemia is the most common etiology for an isolated abducens nerve palsy, but in certain cases neuroimaging is warranted to prevent life-threatening complications. This case highlights the importance and urgency of identifying and managing such conditions.

Sparse Deep Neural Network for Encoding and Decoding the Structural Connectome.

Brain state classification by applying deep learning techniques on neuroimaging data has become a recent topic of research. However, unlike domains where the data is low dimensional or there are large number of available training samples, neuroimaging data is high dimensional and has few training samples. To tackle these issues, we present a sparse feedforward deep neural architecture for encoding and decoding the structural connectome of the human brain. We use a sparsely connected element-wise multiplication as the first hidden layer and a fixed transform layer as the output layer. The number of trainable parameters and the training time is significantly reduced compared to feedforward networks. We demonstrate superior performance of this architecture in encoding the structural connectome implicated in Alzheimer's disease (AD) and Parkinson's disease (PD) from DTI brain scans. For decoding, we propose recursive feature elimination (RFE) algorithm based on DeepLIFT, layer-wise relevance propagation (LRP), and Integrated Gradients (IG) algorithms to remove irrelevant features and thereby identify key biomarkers associated with AD and PD. We show that the proposed architecture reduces 45.1% and 47.1% of the trainable parameters compared to a feedforward DNN with an increase in accuracy by 2.6 % and 3.1% for cognitively normal (CN) vs AD and CN vs PD classification, respectively. We also show that the proposed RFE method leads to a further increase in accuracy by 2.1% and 4% for CN vs AD and CN vs PD classification, while removing approximately 90% to 95% irrelevant features. Furthermore, we argue that the biomarkers (i.e., key brain regions and connections) identified are consistent with previous literature. We show that relevancy score-based methods can yield high discriminative power and are suitable for brain decoding. We also show that the proposed approach led to a reduction in the number of trainable network parameters, an increase in classification accuracy, and a detection of brain connections and regions that were consistent with earlier studies.

Early Neonatal Presentation and Neuroimaging of Parechovirus Meningoencephalitis in a Preterm Baby: A Case Report.

Neonatal meningoencephalitis caused by human parechovirus infection is being increasingly recognized in recent literature. While most cases are postnatally acquired, intrauterine infection is rare, presents early and has a more severe impact on brain health and development. We discuss here an infant born preterm at 34 weeks gestational age, with neonatal course remarkable for severe encephalopathy presenting on day 2 of life due to human parechovirus meningoencephalitis transmitted in utero. Early magnetic resonance brain imaging detected extensive white matter injury and subsequently evolved into multicystic leukoencephalopathy. Posthospital discharge, infant was noted to have early neurodevelopmental impairment at 4 months corrected age.

A novel biomarker selection method using multimodal neuroimaging data.

Identifying biomarkers is essential to obtain the optimal therapeutic benefit while treating patients with late-life depression (LLD). We compare LLD patients with healthy controls (HC) using resting-state functional magnetic resonance and diffusion tensor imaging data to identify neuroimaging biomarkers that may be potentially associated with the underlying pathophysiology of LLD. We implement a Bayesian multimodal local false discovery rate approach for functional connectivity, borrowing strength from structural connectivity to identify disrupted functional connectivity of LLD compared to HC. In the Bayesian framework, we develop an algorithm to control the overall false discovery rate of our findings. We compare our findings with the literature and show that our approach can better detect some regions never discovered before for LLD patients. The Hub of our discovery related to various neurobehavioral disorders can be used to develop behavioral interventions to treat LLD patients who do not respond to antidepressants.

Spontaneous Intracranial Hypotension: Clinical Presentation, Diagnosis, and Treatment Strategies.

Spontaneous intracranial hypotension (SIH) typically presents as an acute orthostatic headache during an upright position, secondary to spinal cerebrospinal fluid leaks. New evidence indicates that a lumbar puncture may not be essential for diagnosing every patient with SIH. Spinal neuroimaging protocols used for diagnosing and localizing spinal cerebrospinal fluid leaks include brain/spinal MRI, computed tomography myelography, digital subtraction myelography, and radionuclide cisternography. Complications of SIH include subdural hematoma, cerebral venous thrombosis, and superficial siderosis. Treatment options encompass conservative management, epidural blood patches, and surgical interventions. The early application of epidural blood patches in all patients with SIH is suggested.

Unsupervised deep representation learning enables phenotype discovery for genetic association studies of brain imaging.

Understanding the genetic architecture of brain structure is challenging, partly due to difficulties in designing robust, non-biased descriptors of brain morphology. Until recently, brain measures for genome-wide association studies (GWAS) consisted of traditionally expert-defined or software-derived image-derived phenotypes (IDPs) that are often based on theoretical preconceptions or computed from limited amounts of data. Here, we present an approach to derive brain imaging phenotypes using unsupervised deep representation learning. We train a 3-D convolutional autoencoder model with reconstruction loss on 6130 UK Biobank (UKBB) participants' T1 or T2-FLAIR (T2) brain MRIs to create a 128-dimensional representation known as Unsupervised Deep learning derived Imaging Phenotypes (UDIPs). GWAS of these UDIPs in held-out UKBB subjects (n = 22,880 discovery and n = 12,359/11,265 replication cohorts for T1/T2) identified 9457 significant SNPs organized into 97 independent genetic loci of which 60 loci were replicated. Twenty-six loci were not reported in earlier T1 and T2 IDP-based UK Biobank GWAS. We developed a perturbation-based decoder interpretation approach to show that these loci are associated with UDIPs mapped to multiple relevant brain regions. Our results established unsupervised deep learning can derive robust, unbiased, heritable, and interpretable brain imaging phenotypes.

Mode-based morphometry: A multiscale approach to mapping human neuroanatomy.

Voxel-based morphometry (VBM) and surface-based morphometry (SBM) are two widely used neuroimaging techniques for investigating brain anatomy. These techniques rely on statistical inferences at individual points (voxels or vertices), clusters of points, or a priori regions-of-interest. They are powerful tools for describing brain anatomy, but offer little insights into the generative processes that shape a particular set of findings. Moreover, they are restricted to a single spatial resolution scale, precluding the opportunity to distinguish anatomical variations that are expressed across multiple scales. Drawing on concepts from classical physics, here we develop an approach, called mode-based morphometry (MBM), that can describe any empirical map of anatomical variations in terms of the fundamental, resonant modes-eigenmodes-of brain anatomy, each tied to a specific spatial scale. Hence, MBM naturally yields a multiscale characterization of the empirical map, affording new opportunities for investigating the spatial frequency content of neuroanatomical variability. Using simulated and empirical data, we show that the validity and reliability of MBM are either comparable or superior to classical vertex-based SBM for capturing differences in cortical thickness maps between two experimental groups. Our approach thus offers a robust, accurate, and informative method for characterizing empirical maps of neuroanatomical variability that can be directly linked to a generative physical process.